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ABSTRACT 

Fluidization velocity is one of the significant parameter to characterize the hydrodynamic studies of fluidized 
bed as it determines different flow regimes. Computational Fluid Dynamics simulations are carried for a cylindrical 
bubbling fluidized bed with a static bed height of lm with 0.150m diameter of gasification chamber. The parameter 
investigated is fluidization velocity in range of 0.05m/s to 0.7mis. Sand with density 2600kg/m 3 and with a constant 
particle diameter of sand 385pm is employed for all the simulations. Simulations are conducted utilizing with the 
commercial Computational Fluid Dynamics software, ANSYS-FLUENT. The bubbling flow regime is appeared above the 
air inlet velocity of 0.2m/s. Bubbling character is increased with increase in inlet air velocities indicated by asymmetrical 
fluctuations of volume fractions in radial directions at different bed heights 

KEYWORDS: Computational Fluid Dynamics (CFD), Bubbling Character, Air Velocity & Volume Fraction 



TRANS 

STELLAR 

•Journal Publications • Research Consultancy 


Received: Sep 12, 2017; Accepted: Nov 02, 2017; Published: Nov 16, 2017; Paper Id.: IJMPERDDEC201738 

NOTATION 

£ s Volume fraction of particle/solid d s Diameter of particle 



Volume fraction of air 

r»max 

Maximum volume fraction of 

particle 

Ps 

Density of particle 

.max 

Maximum viscosity of particle 

p g 

Density of air 

c D 

Standard drag coefficient 

Vs 

Velocity of solid (vector) 

Re s 

Particle Reynolds number 

t s 

Stress tensor for solid 

u mf 

minimum fluidization velocity 

tg 

Stress tensor for gas 

uf 

Advection velocity of solid 

4 

Bulk viscosity of gas 

< 

Advection velocity of gas 

Ps 

Solid pressure 

hs 

Drag force 
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Pf 

Frictional pressure 

Ms 

Solid phase viscosity 

Pk 

Kinetic pressure 

M f 

Frictional viscosity 

Pc 

Collisional pressure 

M k 

Kinetic viscosity 

P g s 

Drag coefficient 

Me 

Collisional viscosity 

So 

Radial distribution function 




e Coefficient of restitution 

INTRODUCTION 

Bubbling fluidized bed is broadly utilized as a part of industrial application due of good blending, warmth and 
mass exchange. Biomass gasification reactor is one for combined heat and power (CHP) production. Flow behavior and 
fluidization properties in the gasifier were studied, by considering the operating parameters such as pressure drop, 
minimum fluidization velocity and bubble behavior. These parameters significantly affect the gasifier efficiency and how 
the hydrodynamic behavior of fluidized bed chamber. Fluidized is critical for deciding on the fine jogging parameters, for 
appropriate fluidization regime [1-4] 

(CFD) was employed to steady the complex phenomena of hydrodynamic actions of gas and solid particles. 
Two models the Eulerian-Lagrangian and Eulerian-Eulerian are accessible in CFD, for the displaying of multiphase 
frameworks. Kinetic theory of granular flow (KTGF) is continuum based completely. This model is greater appropriate for 
simulating big and complicated commercial fluidized mattress. Fluidized bed chambers containing, billions of solid 
particles. In precept, discrete particle models (DPM) is deliver such facts [5-7]. For exceedingly small gas feed the 
chamber may additionally contain a dense bed of fluidized solid particles. The matters can be homogeneously fluidized or 
gas may pass via the bed in the shape of huge bubbles. Further increasing in gas flow rate decreases the bed density and the 
gasoline - solid containing sample may alternate from dense mattress to turbulent mattress then to fast-fluidized mode and 
in the end to pneumatic conveying mode. Those fashions however require cords about strong phase rheology and particle - 
particle interaction legal guidelines. In most of these float regimes the relative importance of fuel - particle, particle- 
particle, and wall interaction is exclusive, it means therefore important to become aware of those regimes to select the 
appropriate mathematical model [5-8]. The vital trouble in modeling of fluidized mattress gasification (FBG) is the 
momentum of levels, that is brief and a massive wide form of independent variables, which include particle density size 
and shape may have an effect on hydrodynamic conduct [2,3,5] 

Taghipour et al. [7], have conveyed exploratory and computational examinations of gas- solid fluidized bed 
chamber hydrodynamics. Their model anticipated bed development and gas-solid steam designs were sensibly well and 
concurring with their experimental results. Huilin et al. [9], has broken down the bubbling fluidized bed with parallel 
blends, applying multi fluid Eularian CFD model, as per the motor hypothesis of granular stream. They found that, the 
hydrodynamics of gas bubbling fluidized bed are related with the distribution of particle sizes and the amount of dissipated 
energy, in particle-particle interactions. Gobin et al.[10], conducted numerical simulation of a fluidized bed, using two- 
phase flow method. They performed time subordinate reproductions, for modern and pilot chamber working conditions. 
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Their forecasts were in great subjective concurrence with the watched conduct as for as bed structure weight drop and 
mean stream administrations. Van Wachem et al.[ll], checked tentatively Eulerian-Eulerian gas-strong model 
reproductions of foaming fluidized beds with existing connections for bubble size or bubble speed utilizing CFX, CFD 
code. The CFD model is based on a two fluid model including the kinetic theory of granular flow. Zhonget. al [12] 
investigated the maximum spoutable bed heights of a spout-fluid bed, filled with six forms of Geldart group D particles. 
The parameters investigated by them are particle size, spout nozzle size and fluidizing gas. They found that, the most 
extreme spoutable bed structure of spout-liquid bed dimensions, with increasing particle size and gas spout size. Fettieri et 
al. [13] simulated the changeover of bubbling to slugging regime for a typical Group B material at four fluidizing 
velocities. The simulations are carried with the Eulerian-Eulerian granular kinetic version available inside the CFX-four 
code. They presented the consequences, in terms of voidage profiles and bubble size, and estimated transition velocity. 
Bahareh Estejab et al. [15], conducted the numerical simulations to predict the hydrodynamic behavior of gas solid 
mixture. They tested different drag models. They found the Huilin-Gidaspow drag model is suitable for single solid phase. 
Embarek Belhadj et al. [16], carried out numerical simulations and experimental validation of hydrodynamics, in bubbling 
fluidized bed combustion. The simulation methodology, adopted by them is Eulerian-Eulerian approach, with KTGF 
theory for solid particles. They validated the results with experimental data and fined to be matched. 

In the present exploration paper, hydrodynamics of two-dimensional non-responsive gas-solid fluidizing 
bedchamber are concentrated computationally. A multi liquid Eulerian model consolidating the kinetic theory for solid is 
applied, in order to simulate the gas-solid stream at various fluidized speeds. 

COMPUTATIONAL METHODOLOGY 

ANSYS FLUENT is utilized for simulation wherein 2 nd segregated first order implicit unsteady solver is utilized 
for multiphase calculations. The usage of Eulerian-Eulerian multi phase version, fashionable k-s dispersed Eulerian multi 
phase model with standard wall functions is used. Gas is taken as continuous phase while sand particles are taken as 
dispersed phase. Kinetic theory of granular flow has been applied to sand particles. Gidaspow model of inter phase 
interaction (Solid-Gas) is used. The bed material considered for the analysis is sand, an inert material. Bed fluidizing 
medium is air. 


The continuity equation for each phase is separately as shown: 



|V s /U+ v 4av)=o 

The only constraint being that total volume fraction has to add up to one. 

£ g +e . = 1 


( 2 ) 


( 3 ) 


In the present work, there is no mass transfer among the phases and thus the terms on the right hand side of the 
equations 1 and 2 are zero. 


The gas phase momentum equation can be expressed as: 
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n)+ v-kp. v ^)= -e,vr + e.P.g + P, k - V,) (4) 

Where P is the pressure, g is the acceleration due to gravity and (3s g is the drag coefficient the stress tensor T g is 
calculated by the following equation: 



Assuming no virtual mass and lift force the solid phase momentum equation can be expressed as: 

^,P,v)+'vk,P,V.v)=V*.-VP,+e.P,g + P ! ,k.-'0 (6) 


Where P is the solid pressure obtained from the Kinetic theory of granular flow? This pressure has three 
components kinetic, collision and frictional. 

SOLVER AND DISCRETIZATION SCHEME 

The section coupled simple technique is performed for pressure-velocity coupling. The second one-order upwind 
scheme is hired for discretization of momentum, turbulence kinetic strength and turbulence dissipation process and the 
primary-order upwind scheme is executed for discretization of quantity-fraction equations. The time step of size is 0.001s 
is taken for the solution to converge. 

GEOMETRY AND MESH 


The geometry parameters of bubbling fluidized bed reactor are primarily based at the studies work of Chinmayee 
Patra [14]. Fig.l (a), shows the geometry of the reactor with its dimensions. The bubbling bed reactor has inner diameters 
of 0.15m and height of lm. The unfastened board region has internal diameters of 0.3m, which shows hydrodynamic study 
and height of 0.8m.The geometry is created by the use of industrial software ICEM CFD. After geometry technology, a 
uniform established mesh has been created, as shown in Fig. 1(b). In this study, total of 39250 cells and 40086 numbers of 
nodes are used, for simulating Fluidized Bed Gasifier. A finer mesh is adapted in the region of static bed and the generated 
mesh, as shown in Fig. (c) 



Figure 1: Schematic of Geometry and Mesh Generated of Bubbling Fluidized Bed 
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RESULTS AND DISCUSSION 


Velocity of fluidizing media, air influence on hydrodynamic behavior of bubbling fluidized bed gasifier has been 
investigated, using CFD simulations. The simulations are carried with CFD commercial software ANSYS FLUENT. 
The results in the form of contour plots, vectors and chats have been presented below for solid volume fraction, phase 
velocity. 

Figure.2 and 3, represents the solid volume fraction contours of sand, for the different air velocities i.e. at 0.05m/s, 
0.2 m/s, 0.5m/s and 0.7 m/s complete and closure view at bed zone. At the velocity of 0.05m/s, the bed begins to expand on 
a small scale, but almost the height of the bed is equal to static bed height. For the air velocity of 0.2m/s, bubbles are 
appeared like surface only and they found in the static bed height, without any appreciable bed growth. The reason 
attributed may be sand debris, within the bottom part of the bed are in pneumatic transport, while fluidization within the 
upper portion is in freely bubbling state. At the air velocity of 0.5m/s and 0.7m/s, it is observed from contour, the bubbles 
flow regime is appeared. This might be reason of the segregate tendencies of the particles, towards the partitions or gulf 
streaming. For that reason, the solid particles slide down along the wall of the reactor, without too much opposition from 
the upward gas flow. 



Figure 2: Sand Volume Fraction Contour at Different Air Velocities for 
Initial Static Bed Height of 0.1m Complete Chamber 



(a) 0.05m/s 


(b) 0.2m/s 


www.tivrc.ors 


editor @tjprc. org 



















346 


B.J.M.Rao, K.V.N.S.Rao & G. Ranga Janardhana 



(c) 0.5m/s (d) 0.7m/s 


Figure 3: Sand Volume Fraction Contour at Different air Velocities for 
Initial Static Bed Height of 0.1m near the Bed Zone 

Flow pattern in fluidized bed is shown in form of the velocity vectors for sand in the Fig.4, for different air 
velocities. It is noticed that, at lower velocities the solid phase velocity is more at the upper surface and lower at the bottom 
i.e., at 0.05m/s as the velocity increased to 0.2m/s, sand is appeared to be lifting in form of surfaces. The bubbling regime 
is come in to sight at the velocity of 0.5m/s. At the velocity 0.7m/s, it is noticed that there is aggressive movement of solid 
particles throughout the bed implying that the velocity at the bottom is less. 
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Figure 4: Velocity Vector of Sand at Different Air Velocities 

Figure.5 and 6 represents the complete and closure views of gas phase velocity vectors in the chamber. The gas 
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phase's velocity vectors in the chamber are heading upwards and velocity is found to be less at the regions of higher solid 
volume fraction as the obstruction is more. For the low velocity 0.05m/s, the air is not able to penetrate properly through 
the sand bed. At the velocity 0.2m/s, it becomes more and increases further with increase in air velocity. 
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Figure 5: Velocity Vector of Air at Different Air Velocities in Complete Chamber 
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Figure 6: Velocity Vector of Air at Different air Velocities in Complete Chamber at Closure View 

Phase 2 i.e. sand volume fractions are plotted against radial length at dissimilar axial heights as shown in Figure.7. 
At air velocity 0.05m/s the bed expansion is at small scale only two heights are considered 0.05m and 0.1m. The volume 
fraction at 0.05m is higher and almost constant in radial direction indicating that low air penetration. In addition small 
scale expansion is indicated at 0.1m bed height with slight fluctuation in sand radial distribution. As the velocity is 
improved to 0.2m/s sand volume fraction at the central core region is less as compared to the wall which indicating the air 
penetration, which indicates higher particle volume fractions along the walls as compared to the core region. With further 
increase in velocities to 0.5m/s and 0.7m/s the bubbling regime appeared as it is indicated by asymmetrical fluctuations of 
volume fractions in radial directions at different bed heights. From the simulation result as shown in the Figure.2. 
The hydrodynamic model is able to describe quantitatively the accumulation of solids near the wall. 
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(c) 0.5m/s (d) 0.7m/s 

Figure 7: Sand Particle Concentration against the Radial Position for 
Different Bed Heights at Different Air Inlet Velocities 

Figure.8 represents the plot of radial variation of axial velocity of sand particles in chamber at different bed 
heights for different air inlet velocities. Sand particle velocities near wall are decreases and this may lead to the back 
accumulation of particles. 




(a) 0.05m/s 


(b) 0.2m/s 


Impact Factor (JCC): 6.8765 


NAAS Rating: 3.11 
































































CFD Analysis of Hydrodynamic Studies of a Bubbling Fluidized Bed 


349 


Ailil «ti»nl at dtimter 




(d) 0.7m/s 


Figure 8: Sand Particle Velocity Distribution along the Radial Direction at Different 
Velocities for Different Axial Heights of 0.05 m, 0.1 m,0.15 m and 0.2m 


CONCLUSIONS 


CFD simulations using ANSYS-Fluent software are carried out to investigate hydrodynamic behavior of bubbling 
fluidized bed. The bubbling flow regime is appeared above the air inlet velocity of 0.2m/s. Bubbling character is increased 
with increase in inlet air velocities indicated by asymmetrical fluctuations of volume fractions in radial directions at 
different bed heights. 
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